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Abstract 

A generalization of the operator product expansion is used to find the dif- 
ferential distributions in the inclusive semileptonic weak decays of heavy 
flavors in QCD. In particular, the double distribution in electron energy and 
invariant mass of the lepton pair is calculated. We are able to calculate the 
distributions in essentially model-independent way as a series in nig 1 where 
niQ is the heavy quark mass. All effects up to nig 2 are included. 
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1 Introduction 



Differential distributions in semileptonic decays of heavy flavors are used for 
measurements of the CKM matrix elements, key phenomenological parame- 
ters of the standard model. To extract the CKM matrix elements from data 
one needs to disentangle the effects of strong interactions at large distances 
from the quark-lepton lagrangian known at short distances. 

Up to now essentially two approaches are applied to describe nonperturba- 
tive strong interaction effects in the inclusive weak decays: the naive parton 
model amended to include the motion of the heavy quark inside the decay- 
ing meson M; and the 'exclusive variant' based on summation of different 
channels, one by one @. Both approaches are admittedly model-dependent, 
neither their accuracy nor the connection to the fundamental parameters of 
QCD are clear a priori. Each of them needs an input from constituent quark 
model to parametrize nonperturbative effects. The latter play an especially 
important role in the form of the spectra near the endpoints. 

The need for the model-independent QCD-based predictions is apparent. 
Considerable progress achieved recently in the theory of preasymptotic effects 
(proportional to powers of I/tuq where rriQ is the heavy quark mass) allows 
one to make these predictions. 

The theoretical construction presented in this paper is, in a sense, a gener- 
alization and combination of the formalisms which are used in deep inelastic 
scattering and total cross section of e + e~ annihilation. The expansion param- 
eter in deep inelastic scattering is Q^ 1 where Q is the momentum transfer. 
In the problem at hand the expansion parameter is nig 1 or, more exactly, 
the inverse energy released in the final hadronic state (in the rest frame of 
the decaying quark). 

In the classical problems of this type, like e + e~ - annihilation, there are 
two alternative ways to get predictions. The first approach having a solid 
theoretical justification in terms of OPE J3J is based on calculations in the 
euclidean domain where one can apply OPE. The contact with the observable 
quantities is made through the dispersion relations and in this way predic- 
tions for certain integrals are obtained. In the second approach we perform 
the calculations directly in Minkowski domain. Although formally this cal- 
culation refers to large distances from the first approach we know that in 
specific integrals large distance contributions drop out. Therefore the results 
obtained in this way, although not valid literally, should be understood in 
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the sense of duality: being smeared over some duality interval the theoretical 
prediction should coincide with the smeared experimental curve. The inclu- 
sive weak decays will be treated within the second approach. The averaging 
mainly refers to the invariant mass of the inclusive hadronic state produced 
in the decay considered. 

If the invariant mass of the final hadronic state is large this is not a 
constraint at all since the theory 'itself takes care of the averaging required 
by duality. In the opposite limit, near a spectral endpoint, the smearing is 
not provided for free. We will discuss the issue in some detail below. 

Although we explicitly work in the Minkowski kinematics we always keep 
in mind the relationship to the euclidean domain and the corresponding op- 
erator product expansion. The first analysis of this type has been outlined in 
[[| for inclusive heavy flavor decay rates. A general analysis of the semilep- 
tonic inclusive spectra along this line is presented in ref. In that work 
it was observed, in particular, that the leading operator and those appear- 
ing at the next-to-leading order have a gap in dimensions of two units, and, 
consequently, the O^uiq 1 ) term should be absent in certain quantities. The 
analysis presented in || was not backed up, however, by concrete calculations 
of the preasymptotic effects. Recently this formalism has been systematically 
developed and applied to the non-leptonic decays of heavy flavors || |6| and 
the charged-lepton energy spectrum in the semileptonic decays |7j (see also 
[§]). The present work is a natural continuation of ref. J?|. 

We generalize the results of ref.]7j to find the complete inclusive distri- 
butions in the semileptonic decays. The leptonic variables — E e ,q 2 and q , 
where E e is the charged lepton energy and q is the momentum of the lepton 
pair f] — are kept fixed which automatically fixes the invariant mass of the 
inclusive hadronic state. Integrating over go we obtain the double spectral 
distribution in E e and q 2 . 

At the first stage we construct the transition operator T(Q —>■ X — > Q) 
describing the forward scattering amplitude of the heavy quark Q on a weak 
current. Our focus is the influence of the 'soft' modes (background fields) 
on the transition operator T^ v which is expressed as an infinite series in the 
local operators built from the gluon and/or light-quark fields bilinear in QQ. 

The local operators are ordered according to their dimensions; the co- 
efficient functions contain the corresponding powers of 1/rriQ (or 1/E, the 

6 The charged lepton produced will be generically called 'electron' hereafter. 
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energy release). At sufficiently large tuq or E the operators with the lowest 
dimensions dominate, and the infinite series can be truncated. 

At the next stage the matrix elements of the relevant operators over the 
initial heavy hadron Hq must be evaluated. Unfortunately, in the present- 
day QCD the matrix elements over the hadronic states are not theoretically 
calculable. In some instances they can be related, through heavy quark 



symmetries, to measurable quantities ||12|| , [g; in other cases they have to 
be parametrized. These parameters play the role analogous to the gluon 
condensate As a matter of fact, at the level of the leading preasymptotic 
corrections only two operators are relevant. The matrix element of the first 
one can be related to the mass splittings of the vector and pseudoscalar 
heavy mesons. The matrix element of the second one has the meaning of the 
average square of the spatial momentum of the heavy quark Q in Hq and 
the state must be treated clS cL parameter. 

Finally, the observed decay rates and spectra are obtained by taking the 
discontinuity of the hadronic tensor < Hq\T^ u \Hq > and convoluting the 
result with the lepton currents and appropriate kinematic factors. 

In this paper we consider the differential distributions in the semileptonic 
decays at the level of 0(rriQ 2 ). The differential distributions are measured 
experimentally in the B meson decays and will be used for more precise 
determination of V u b, for example. This was a primary motivation for our 
investigation. We would like to make it as close to the fundamental QCD as 
possible. 

The organization of the paper is as follows. In Sect. 2 we describe the 
kinematics and in Sect. 3 we present the operator product expansion. In 
Sect. 4 we derive the differential distributions. Sect. 5 is devoted to the 
analysis of our distributions and limitations on the their use. Our results are 
summarized in Sect. 6. Appendix contains expressions for hadronic invariant 
functions. 



2 Kinematical analysis 

We will consider the inclusive weak decays of the mesons (or baryons) with 
the open heavy flavor into the lepton pair plus (inclusive) hadronic state 
Hq(ph) l(pi) + v(p u ) + hadrons. Our final goal is to calculate the differ- 
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ential decay rate 

d 3 T 

— (1) 

dE e dq 2 dqo 

where E e is the energy of the emitted electron and g M = pf + p^ is the 4- 
momentum of the lepton pair. In order to find the differential distributions we 
need to know the amplitude of the process, which is given by the expression 

M = V qQ ^eV v v <X\j v \H Q > (2) 

Here V q Q is the corresponding Cabibbo-Kobayashi-Maskawa matrix element, 
j M = qT^Q is the electroweak currents, T M = 7 M (1 + 75). (Although our 
theory is general we will keep in mind the b — *> c and b — > u decays, so that 
Q = b and q = c or u). The differential distributions we are interested in are 
given by the modulus squared of the amplitude (0) summed over the final 
hadronic states. 

The modulus squared of the amplitude summed over the final hadronic 
states can be written as 

\M\ 2 = \V qQ \ 2 G 2 F M Hq l^W, u (3) 

where Mh q is the mass of hadron Hq , W^ v is the hadronic tensor 

WV = W5\p Hq -q-p x )Y: — ^- < H q \jI(0)\X X X\U0)\H Q >, 

(4) 

and V" is the lepton tensor 

l>» = 8 [ {p e Y{p v ) V + (PeYiPuT ~ <T (Pe " Pv) + (5) 

Let us introduce the hadronic structure functions Wi and parametrize the 
hadronic tensor in the following way: 

Wp„ = -Wi g^ v + w 2 v^v u -i w 3 e^ a p v a qp + 

Wi q^ q v + w 5 {q u v^ + q^v u ). (6) 

Here q^ = (p e +Pu)fi is the 4 - momentum of the lepton pair, t> M = {j> Hq )^/ M Hq 
is the 4-velocity of the initial hadron. Note that we have omitted the structure 
q^Vy — qyV^ which can not appear because of the T-invariance. The structure 
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functions Wi depend on two invariant variables, q-v and q 2 . In the rest frame 
of Hq which will be used throughout the paper q ■ v = q , so Wi = Wi(qo, q 2 ). 
The convolution of W^ v with the lepton tensor (|5|) is given by the expression: 

W^l^ = 4{2q 2 w 1 + [4E e (q -E e )-q 2 }w 2 + 2q 2 (2E e -q )w 3 }. (7) 

We see that only three invariant functions are relevant for the processes we 
are considering in this paper. At this step we encounter the third variable, 
the electron energy E e = p e - p HQ /M Hq , entering through the leptonic tensor. 
Finally the formulae for the differential width takes the form 

d 3 T G 2 

= [2<A °' + [4 E ' Xq ° ~ ^ ~ " 2] W2 + 2 " 2(2 E ' ~ *> "> 1 

(8) 

This expression concludes the kinematical analysis. Our task is, of course, 
the calculation of the invariant functions Wi(q ,q 2 ). We will proceed to this 
calculation in the next section. 



3 Operator product expansion 

In this section we will discuss the derivation of the tensor W^. The Opera- 
tor product expansion is similar to that in the deep inelastic scattering. It is 
convenient to introduce the hadronic tensor (forward scattering ampli- 
tude) 

h lv = ij dw^^- < H Q \T{j;{x)j v m\H Q > . (9) 

The absorptive part of this tensor reduces to discussed above 

WV = (l/z)disc(V)- (10) 

Here disc(/i M ^) is the discontinuity of the forward scattering amplitude h^ v 
on the physical cut in the complex plane of the variable go- Of course, h^ v 
can be expanded into the same set of structures as W^ v (see. eq. (^j)) 

h^ u = -h g^ v +h 2 Vfj,v u -i h 3 e^ a/3 v a q p + h 4 q^ q v + h 5 (q v + g M v v ), (11) 

and the relation (|H]) implies that 

w l = 2\m{h i ). (12) 
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Let us remind that h^ u is the matrix element of the transition operator 

T 

T, v = i J d 4 xe~^ T{j+(x)j v (Q)}, (14) 

so below we will construct OPE for the product of currents in equation fll^p. 
Having in mind the relationship to euclidean analysis discussed above we 
will treat our expansion in the same way as a normal euclidean OPE. In the 
asymptotic limit itlq — > oo the hadronic tensor is given by the tree graph 
of Fig. 1. This graph defines the matrix element of the transition operator 
T^y over the heavy quark state, 

< Q\Tfu,\Q >= -uqT,—^- T v u Q (15) 

if— q -m q 

This expression represents nothing else but the free quark decay. In the 
asymptotic regime tuq — > oo the interaction of the heavy quark with the 
gluon/light quark medium, as well as its intrinsic notion inside the hadron 
can be neglected. Then 

Pfj. = Pov = rriQVf, (16) 

where is the Hq 4-velocity. 

Equation fll~5|) allows one to immediately write down the operator form in 
the approximation at hand (only operators bilinear in Q, Q are considered, 
see a discussion of other operators at the end of the section): 

P^u = -QT ft -^—T v Q = (17) 



- m q 



(k 2 — mT) 



where k = Pq — q. 

As we see, the two operators Q^ a Q and Q^albQ showed up in the operator 
expansion at the level considered. Note, that the Q^alsQ term vanishes after 
averaging over the unpolarized hadronic states. 

In this paper the perturbative corrections in a s are not touched upon at 
all. As for nonperturbative corrections they appear due to interactions with 
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the soft medium of the light cloud in Hq. By taking these interactions into 
account we isolate two types of effects. First, the fast quark q produced does 
not propagate as a free one, but interacts with the background fields; these 
corrections will be included explicitly. Second, the heavy quark Q also does 
not live in the empty space; it is surrounded by the light cloud. In particular, 
due to this fact the heavy quark momentum does not coincide with rnQV^. 
This large distance effect will not be calculated explicitly, but implicitly it 
will be reflected in the Hq matrix elements of the operators in T^ v . This is in 
full analogy with what people usually do in the deep inelastic scattering. The 
influence of the background field on the transition operator is summarized 
by the following expression 

T ^ = ~J dxe- iqx Q(x)T^S q (x, 0) Y v Q(0) (18) 

where S q (x, 0) is the propagator of the quark q in an external gluon field A a . 
It is convenient to use the Schwinger technique of treating the motion in an 
external field (for a review of QCD adaptation see , e.g. ref. [0). Within 
that formalism the propagator S q is presented by the following expression 

^(a., 0) = -1— -|0) (19) 

If Tflq 

Here J> = 7^(p M + A^(X)), = g A a ^T a is a gluon field in the matrix 
representation. Furthermore, the operator of coordinate and momentum 
Pn are introduced, (thus the field A^(X) becomes an operator function of 
Xfj), with the commutation relations 

[ P ^X U ] =ig^ [X^X V \ = 0, \p^p v \ = 0. (20) 

The states \x) are the eigenstates of the operator X^, X^\x) = x^\x). 
Combining equations (|T9|) and (|1SD , we arrive at 

T IUI = - [ dxe-^(x|Q(x)r M — ^- T u Q(X)\0). (21) 

As we have discussed above the operator contains a large mechanical 
part (-Po)^t = m QV^, the deviation from P will be separated explicitly 

V» = (PoV + ^ (22) 
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and we will expand in ir^. In this paper we will limit ourselves to the terms 
up to 0(tt 2 ) corresponding to l/mn corrections. The master formulae to 
perform the expansion is 

T^ = - f dx(x\Q(X)T,- -± —r„Q(X)\0). (23) 

There is a subtle point in the description of the formalism given above. 
Technically in the computation the A^{x) is assumed to be a c- number back- 
ground field while in the final expression for local operators it should be 
understood as a second quantized operator. Since we are not considering 
any loop corrections this substitution is justified. 

Let us now discuss the set of the operators relevant to the order 0(rriQ 2 ). 
Without loss of generality we can work in the rest frame of the hadron Hq, 
i.e. = (1,0,0,0). Only those operators will be retained which produce 
non- vanishing results after being averaged over Hq. The leading operator, 
as it was discussed above, is 

QloQ, (24) 
it's matrix element is fixed by the vector current conservation to be 

— l — < H Q \Q l0 Q\H Q >= 1, (25) 
1 m h q 

Equation ( fffi) is given in relativistic normalization we are using throughout 
this paper. In the non-relativistic normalization there is no need in the factor 
1/2M Hq in the l.h.s. 

As it has been noted in the ref. || there are no operators of dimension 4 
in the problem at hand. The set includes two operators of dimension 5: 

O g = ^Qcr aP G a/3 Q, (26) 

O n = -QD 2 Q = Qn 2 Q, (27) 

where a a)3 = ~ (7°7 /3 — 7 /3 7 Q ), and G a p = g G q a/3 T a is the gluon field strength 
tensor. The classification above takes into account the fact that the quark 
field Q satisfies the equation of motion. In particular, it stems that the 
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operator QQ is not independent but is reducible to three operators ([24]) , 
(|D and (|7|): 

QQ = QloQ- 7 ^ T Q^Q + J 1 T Q^ G ^Q + 0(m Q 3 ). (28) 
Z vtlq 4 mg 

To get eq. (28) we observe that the lower component of Q is related to the 
upper one in the following way 

— Q = ^r™— Q + ° (m « } ' (29) 

and the difference between QQ and Q-yoQ is due to the product of the lower 
components. (Here and below we will stick to the Hq rest frame.) 

A few other useful relations which can be obtained in the same manner 
and are valid at the level 0(rriQ 2 ) are: 

QinQ = -^Q(7f 2 - t -aG)Q + 0(mQ 2 ), (30) 

Q^ lo Q = 0(mQ 2 ), (31) 
Q*oQ = ^Q(* 2 - \ aG )Q + 0{m Q 2 ). (32) 

A few comments are in order here concerning the actual technique of 
constructing the OPE. Since we work in the Hq rest frame it is convenient 
to compute different components of T M „ separately, T 00 , T oi , T i0 and T^. 
The calculation itself is a straightforward although rather tedious procedure 
of expanding the denominator in eq. (^) in jt using the properties of the 7 
matrices, the commutation relation 



i (33) 



and equations ( |30"D - (|32j). 

Notice that we must keep the terms of the first order in tt and of the 
second order in 7?, since 

n Q= { -^Q + O(m Q 2 ). (34) 
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Next, observe that the Green function in the background field can be written 
as follows: 



f-4-m q w qj (V-q) 2 + (t/2)aG-m 2 q w 9/ lT 

(35) 

To transpose 1/11 with T u it is convenient to use the identity 

Ir„ = rl + I[r„,n]i 
= r„I + I|T„i ff GlI. (36) 

Acting on Q and using the equations of motion we can now substitute 1/11 
in both terms on the right-hand side by 

1 

m 2 Q -m 2 q - 2Vq + q 2 ' ^ 

provided that we limit ourselves to terms up to 0(mg 2 ). The second term 
in ( |3~6|) can be simplified even further since here we can additionally neglect 
Ti in V = Po + 7T. 

We split the calculation into three parts: Vector x Vector, Axial x Axial 
and Vector x Axial in correspondence with the structure of as a sum of 
vector and axial vector, = 7^ + 7^75. The full hadronic tensor h^ u is given 
then by the following expression: 

h — h VV J-h AA J-h AV J-h VA — 1 / IT \rpVV .rpAA.rpAV j_rpVA\ TT ^ 

- n fiu + n ^u + tl f i» + tl l iv - 7^ < n Q\ I f iu + I f iu + 1 t xv \ n Q > ■ 

(38) 

The complete expressions for the hadronic invariant functions are given in 
the Appendix. In the order 0(rriQ 2 ) they are defined by matrix elements of 
operators Oq, O n given by eqs. ([26|) , (|7|): 

' <H Q \Q t -a IMU G^Q\H Q >= f i 2 G , (39) 



2M Hq 2 



' <H Q \Q7f 2 Q\H Q >=^, (40) 



10 



The parameter y? G coincides with m 2 H introduced in |TJ]]. For mesonic states 
it is expressible in terms of the quantity measured experimentally - the hy- 
perfine mass splittings, and it has the zero value for baryonic states of the 
type of Aq; fi 2 has the meaning of the average square of spatial momentum 
of the heavy quark Q in the hadronic state Hq. The two parameters, [Xq and 



^1, often appear in the combination fi 2 — /i^., cf. equation (p8|). 

The last comment of this Section is about the operators which are not 
bilinear in Q, Q fields. The simplest example of appearance of such operators 
is given by the diagram of Fig. 2 where the heavy quark Q propagates between 
the current vertices. This diagram is similar to the one of Fig.l, and the 
corresponding operator follows from eq.(pT8|) by substitution Q, Q => q, q, 
m q =^ tuq, =^> q^. The additional term in has the form 

AT^ = -qT v — - T^q= (41) 

2 



(q 2 - m\ 

Matrix element of the operator (pi a q over Hq state counts the number 
of quarks q and is not small in general. The operator coefficient given by 
eq.(|42|) is particularly large when q 2 — > tuq. 

In terms of the intermediate hadronic states in the forward scattering off 
Hq this contribution is due to states in crossing channel containing two Q 
quarks - the corresponding problem was pointed out in ref. 0. The cross- 
channel is not related to the weak inclusive decays under consideration and 
the analysis above shows that we can consistently omit the crossing channel 
together with operators in related to this channel. Similar consideration 
can be carried out for other operators containing no heavy quark fields Q, Q. 



4 Calculation of the differential distributions 

The differential distributions we are interested in are determined by equation 
(^) containing three invariant functions w±, w-i and W3. They are obtained 
from the results for hi (see Appendix) by taking the imaginary parts of the 
corresponding functions (see. eq. flT2|)). The imaginary parts are due to the 
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poles of hi and are obtained through the relations: 



1 (—l)n-l An-l 

Im- = 7r l) ^-—6(z) (42) 



where z is given by 



2 Q -2m Q q + q 2 - mj. (43) 



We don't present here the expression for the triple differential distribution 
which can be easily obtained by combining the equations (§), (|I2|) , (f42"D and 
expressions for hi from the Appendix. 

Although the result is written for the physical quantity d 3 T/dE e dq 2 dqo, 
it can not be directly compared with the experimental data. An obvious 
signal for this is the presence of the delta function and it's derivatives. It 
is not surprising because we are sitting now right on the physical cut on 
mass-shell of the g-quark. As we discussed in the introduction our results 
should be understood in the sense of duality: that is that the predictions 
should be smeared over certain duality interval. At the moment we have 
no purely theoretical tools to fix the size of the duality interval, therefore 
we are forced to rely on qualitative arguments and experimental data. For 
example the duality interval for q Q can be inferred from the distribution in 
the invariant mass of the final hadronic states. Our (^-functions reflect the 
resonance structure at low invariant masses. The smearing interval should 
be chosen in such a way as to cover the entire resonance domain up to the 
onset of the smooth behavior. Instead of smearing of the distribution one can 
calculate the average characteristics like the total width T or < M x >, where 
Mx is the invariant mass of the final hadronic states. The power corrections 
we have calculated will enter in a specific way in each particular quantity. 

Now let us proceed to the calculation of the double differential distribution 
d 2 T / dE e dq 2 . To this end we must integrate over go, rather simple exercise 
with (^-functions. However if one would perform the integration by merely 
substitution 

m% + q 2 — m 2 
\m Q (44) 

and taking the derivatives in the case of 5' and 5", one would get the wrong 
answer. The point is that integration domain in q has a boundary from 
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below 

qo>E e + ^L_, (45) 

which corresponds to 4E e E u > g 2 . Therefore one should take into account 
the fact that go can n °t cross the boundary p5[) . For that we introduce 
#(go — E e — q 2 /AE e ) into the integrand. The occurrence of the ^-function is 
important for the integration of 5'(q — q£) and S"(q — q^) which leads to 
appearance of 5(qQ — E e — q 2 /4E e ) and 8'{ql — E e — q 2 /4E e ) in the double 
distribution d 2 T /dq 2 dE e , because of differentiation of the function . The 
final formulae for the double differential distribution in the lepton energy E e 
and q 2 takes the form: 

G Q [ 1 - 5p + 2t + 10pt + lOart - 10xt 2 - 
(-1 + Qp - hp 2 + x - 5px + t - 2pt + 5p 2 t + xt + 15pxt + 
hxH - 2xt 2 - lOpxt 2 - 10x 2 t 2 + 5ar¥) -t)(l-x)-p)] + 

K Q [-3 + 3p + 4t - 4pt - 6xt + Axt 2 - 
(1 - 2p + p 2 - 3x + 3px -3t + 2pt + p 2 t + llxt - 3pxt - 
3x 2 t - Qxt 2 - 2pxt 2 + 2x 2 t 2 + x 2 t 3 ) - t)(l - x) - p) + 
(1 - p - x + xt){l - t){\ - 2p + 

p 2 - 2xt - 2pxt + x 2 t 2 ) -t){l-x)-p)\ }. (46) 
Here we have introduced the dimensionless variables 

x = 2 E e /m Q , t = q 2 /2m Q E ej (47) 

and the parameters 

p = m 2 q /m 2 Q , G Q = p 2 G /m 2 Q , K Q =pl/m 2 Q . (48) 

Let us emphasize that the scale ttiq used in equation (|7|) is the heavy quark 
mass and does not coincide with Mh q which is normally used in the experi- 
mental distributions. 

The fact that OPE generates corrections only of the order of 0(mq 2 ) 
(terms proportional to Kq and Gq) is valid for the distributions only if we 
use m,Q as a scale, i.e. in the variables x, t. Of course one can easily rescale 
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them to Mh q ] then the corrections of the order of 0(rriQ ) will show up for 
trivial kinematical reasons. 

We can proceed further and obtain the energy spectrum by integrating 
over q 2 . The range of integration is given by 

< t < 1 — . (49) 

1 — x 

The result for the energy spectrum coincides with that obtained in 0. For 
the sake of completeness we present it here []: 

fx = |VW|2 T9^ (1 ~ X ~ P)2x2{(1 " /)2(1 + 2/)(2 ~ X) + (1 " /)3(1 ~ X) + 

Pi If) f 2 

(1 - /)[(1 - /)(2 + -x-2f + —fx) - J -(2x + /(12 - 12x + 5x 2 ))]G Q - 

[ 1(1 - f) 2 (l + 2f)x + £- (1 - f)(10x - 8x 2 ) + 

6 p 

ll(3-4f)(2x 2 -x 3 )}K Q }, (50) 
P 

where f — p/(l — x). Finally, performing the last integration over x in the 
domain 

0<a;<l-p, (51) 
we arrive to the total width coinciding with that in ||: 

r = l^| 2 ^?[^o (1 + \ (G Q - K Q )) -2 Zl G Q ], (52) 

where z = 1 — 8p + 8p 3 — p 4 — 12p 2 log p and Zi — (1 — p) 4 . 

Now let us discuss the characteristic features of the double distribution 
(PI). The most striking one is the presence of the singular terms. The 
technical reason for occurrence of those terms was that we expanded the 
denominator of the pole expression (|2"3| ) in n and a G. Physically this expan- 
sion reflects the shifts of the masses of particles due to the nonperturbative 
effects. As it was mentioned above these singularities reflect the structure 
of the resonance domain and the predictions suitable for comparison with 



7 Let us draw the reader's attention to the difference of notation: y in Q is equal to 
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the experimental data require smearing over corresponding interval. To il- 
lustrate the most salient features of our prediction let us concentrate on the 
physically interesting case of the b — > u transition. 

For massless u quark the kinematical region of h quark semileptonic decay 
is shown on Fig. 3. It has the form of a square with the side equal to 1 
in the plane (x = 2E e /mb, t = q 2 /2rribE e ). The right-hand side of the 
square corresponds to the maximal energy of electron E e = rrib/2 while the 
upper side is a maximal energy of neutrino. In the real B meson decay 
the kinematical region is certainly wider; if one neglects the pion mass the 
region is the square with the side x max = t max = Ms/rrib. The origin of this 
window is related to the motion of the heavy b quark inside of B meson. 
In our calculations we account for nonzero momentum of the b quark in the 
form of expansion which produced singular 5 and 8' terms on the boundary. 
It is possible to show (see refs.P], that the expansion breaks down at 



distances ~ (M B — m b ) / m b near the boundary, so we need to average our 
results over the range of the order of the window between quark and hadron 
boundaries. It is interesting to note that the distribution spread off the 
distances of the order (Mb — rrib)/mb while the corrections to integrals are 
only of the second order in l/m&. 

Another effect we need to account for is the structure of resonance region 
near the low end of the hadronic invariant masses. To imitate the effect let 
us imagine that this region corresponds to the u quark fragmentation into 
the hadronic states with s (the square of the invariant mass) from s = to 
s = s = 2GeV 2 . The curve corresponding to s = s Q at Fig. 3 is given by the 
equation: 

(l-t)(l-x) = s /m 2 b (53) 

and the resonance region should be included as a whole into the process of 
averaging; we can predict the integral but not the structure. 



5 Application to the analysis of the experi- 
mental data 

In order to compare our results with the experimental data we use the fol- 
lowing values of parameters in eq. (^6j). First, we use m b ~ 4.8 GeV as 
deduced from QCD analyses of the Ypsilon system and M B ~ 5.3 GeV. For 
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the parameter Gb we use the value [|l4f 



G b = -; 




(54) 



As a representative value we use for the parameter Kb the value ~ 0.02 (see 



To construct the quantity convenient for the comparison of the theory 
with experiment we consider the integral 



where x c , t c is the point in (x, i) plane sitting not too close to the boundary 
(outside the resonance range), T = \V u b\ 2 G 2 F m\/ 1^2^ and A(x c ,t c ) is the 
area of integration shown on Fig. 3 which includes the resonance domain. 

The function P(x, t) is plotted as a function of t on Fig.4 for three values 
of x equal to 0.3, 0.6, 0.8. The last value of x is close to the border of the 
resonance region beyond which we cannot make reliable predictions for the 
distributions considered. The dashed lines on Fig.4 describe the leading order 
distributions in t while the solid lines include QCD corrections we calculated. 
As we can see it from the curves, the corrections are negative. 

6 Conclusion 

Let us now summarize our results. Model independent approach to nonper- 
turbative effects {l/m q ) n is used for calculations of differential distributions. 
The effects are most pronounced near the endpoints of the spectra. We dis- 
cussed how the comparison with experiment should be formulated accounting 
for the boundary effects. A somewhat disappointing point is that we cannot 
use our results to improve an extraction of V u b by the consideration of q 2 
dependence. Indeed, experimentally the signal of b — > u is due to the range 
of electron energy E e near the upper end where b — ► c is absent. However as 
it follows from Fig. 3 the distribution in q 2 at such energies is concentrated 
in the resonance domain, and no model-independent prediction emerges. 



ref. H-B). 




(55) 
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A Hadronic invariant functions 



Here we present the results of calculations of different hadronic invariant 
functions hi, introduced by eqs. flTTl , (|38D ). The structure functions Wi are 
simply related to hi by eqs. (|12"D and (f42|). We use the following notation : 
go = q • v, q 2 = q 2 — q 2 and z = ttiq — 2 ttiq q + q 2 — m 2 . 
For the Vector x Vector functions we have: 



— [lv 2 G (( 4m Q-3q )(m Q -m q -q ) + 2q 2 ) + 

m Q "J 



A (<?o (m Q - m q - q ) - - q 2 > i — 



2 .2X1 1 

3 9 )] ^ 



^ /4 g 2 (m Q - m 9 - g ) ^ , (A.l) 



2 i 

- [2// G {m Q -m q ) - 5/i G g + 7/4<?o] 7^ 



^ 2 ^2 1 

^rn Q ^q - 



mgfilq 2 ^ (A.2) 



^ = 0, (A.3) 



^r = -^(/4-/4)-, (a.4) 

oTfln Z 



!S ---^[5^ G -^)-4^)]i + |^g 2 i 

2 3 TOg 3 Z^ 



To get the functions hf A for Axial x Axial tensor from hj v one should 
substitute m q by (— m g ) in eqs.( |A.l| - |A.5|) . 

For the Axial x Vector tensor only one invariant structure survives: 



z 3 niQ z l 3 z 6 
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Summing up we get the result for the full hadronic tensor h^ v . 

h = -[2(m Q -q ) - 7^- {Hg-&)\\- 

Q 

[^/4 (4m 2 Q + 2q 2 - 7 m Q q + 3q 2 ) + 
" (4go(m Q -g )-^ 2 )]^- 



Ivtiq 3 z 1 

!| /4 9 2 (™Q ~ <?o) , ( A -7) 



10 / 2 2m 1 

h2 = -[^m Q + — ^-^ a )]-- 



.28 2 2 .8 20 . . 1 

[— ^q + fi G (-m Q - —q )\^ - 

— £m Q q 2 -. (A.8) 

^--db"*-"^ < Aio) 

fc5 = 2i-|[5(^-ri)^-4ri)]i + |ri f 2 i (A.11) 
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Figure captions 

Fig.l: The tree diagram determining the transition operator T^ v in the 
leading approximation. The dashed lines correspond to the weak currents, 
the solid internal line describes the propagation of the quark q and the bold 
external lines represent the heavy quark Q. 

Fig. 2: The tree diagram determining the operator without the heavy 
quark Q. Now the bold internal line describes the propagation of the heavy 
quark Q and the solid external lines represent the quark q. 

Fig. 3: The kinematical region of the decay for b — > u decays in coordi- 
nates x = 2E e /mb and t = q 2 /2m b E e . The solid lines are the kinematical 
boundary for the b quark decay (x max = t max = 1) and the dashed lines are 
the boundary for B meson decay (x max = t max = Mb /nib). The area of 
integration for the distribution P(x,t) (see eq.flsB])) is shadowed. It includes 
integration over the resonance domain. 

Fig. 4: The integrated distribution P(x, t) (see eq. (|55"D ) for the case b — > u 
is plotted as function of t = q 2 /2m b E e for few values of x = 2E e /m b . The 
dashed lines correspond to the leading order distribution while the solid lines 
account for QCD corrections. The lines stop at the border of resonance 
region. The corrections are negative and their relative magnitude becomes 
larger when x approaches the resonance region. 
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